SaLTy: a novel Staphylococcus aureus Lineage Typer.

Staphylococcus aureus asymptomatically colonises 30 % of humans but can also cause a range of diseases, which can be fatal. In 2017 S. aureus was associated with 20 000 deaths in the USA alone. Dividing S. aureus isolates into smaller sub-groups can reveal the emergence of distinct sub-populations with varying potential to cause infections. Despite multiple molecular typing methods categorising such sub-groups, they do not take full advantage of S. aureus genome sequences when describing the fundamental population structure of the species. In this study, we developed Staphylococcus aureus Lineage Typing (SaLTy), which rapidly divides the species into 61 phylogenetically congruent lineages. Alleles of three core genes were identified that uniquely define the 61 lineages and were used for SaLTy typing. SaLTy was validated on 5000 genomes and 99.12 % (4956/5000) of isolates were assigned the correct lineage. We compared SaLTy lineages to previously calculated clonal complexes (CCs) from BIGSdb (n=21 173). SALTy improves on CCs by grouping isolates congruently with phylogenetic structure. SaLTy lineages were further used to describe the carriage of Staphylococcal chromosomal cassette containing mecA (SCCmec) which is carried by methicillin-resistant S. aureus (MRSA). Most lineages had isolates lacking SCCmec and the four largest lineages varied in SCCmec over time. Classifying isolates into SaLTy lineages, which were further SCCmec typed, allowed SaLTy to describe high-level MRSA epidemiology. We provide SaLTy as a simple typing method that defines phylogenetic lineages (https://github.com/LanLab/SaLTy). SaLTy is highly accurate and can quickly analyse large amounts of S. aureus genome data. SaLTy will aid the characterisation of S. aureus populations and ongoing surveillance of sub-groups that threaten human health.

A Multilocus Sequence Tying Scheme for Rapid Identification of Xanthomonas citri Based on Whole Genome Sequencing Data.

Xanthomonas citri is a plant-pathogenic bacterium associated with a diverse range of plant host species. It has undergone substantial reclassification and currently consists of fourteen different subspecies or pathovars that are responsible for a wide range of plant diseases. Whole genome sequencing (WGS) provides a cutting-edge advantage over other diagnostic techniques in epidemiological and evolutionary studies of X. citri because it has a higher discriminatory power and is replicable across laboratories. Also, WGS allows the improvement of multi-locus sequence typing (MLST) schemes. In this study, we used genome sequences of Xanthomonas isolates from the NCBI RefSeq database to develop a seven-gene MLST scheme that yielded 19 sequence types (STs) that correlated with phylogenetic clades of X. citri subspecies/pathovars. Using this MLST scheme, we examined 2,911 assemblies from NCBI GenBank and identified 15 novel STs from 37 isolates that were misclassified in the NCBI. In total, we identified 545 X. citri assemblies from GenBank with 95% average nucleotide identity to the X. citri type strain and all were classified as one of the 34 STs. All MLST classifications correlated with phylogenetic position inferred from alignments using 92 conserved genes. We observed several instances where strains from different pathovars formed closely related monophyletic clades and shared the same ST, indicating that further investigation of the validity of these pathovars is required. Our MLST scheme described here is a robust tool for rapid classification of X. citri pathovars using WGS and a powerful method for further comprehensive taxonomic revision of X. citri pathovars.

Uncovering the boundaries of Campylobacter species through large-scale phylogenetic and nucleotide identity analyses.

Campylobacter species are typically helical shaped, Gram-negative, and non-spore-forming bacteria. Species in this genus include established foodborne and animal pathogens as well as emerging pathogens. The accumulation of genomic data from the Campylobacter genus has increased exponentially in recent years, accompanied by the discovery of putative new species. At present, the lack of a standardized species boundary complicates distinguishing established and novel species. We defined the Campylobacter genus core genome (500 loci) using publicly available Campylobacter complete genomes (n = 498) and constructed a core genome phylogeny using 2,193 publicly available Campylobacter genomes to examine inter-species diversity and species boundaries. Utilizing 8,440 Campylobacter genomes representing 33 species and 8 subspecies, we found species delineation based on an average nucleotide identity (ANI) cutoff of 94.2% is consistent with the core genome phylogeny. We identified 60 ANI genomic species that delineated Campylobacter species in concordance with previous comparative genetic studies. All pairwise ANI genomic species pairs had in silico DNA-DNA hybridization scores of less than 70%, supporting their delineation as separate species. We provide the tool Campylobacter Genomic Species typer (CampyGStyper) that assigns ANI genomic species to query genomes based on ANI similarities to medoid genomes from each ANI genomic species with an accuracy of 99.96%. The ANI genomic species definitions proposed here allow consistent species definition in the Campylobacter genus and will facilitate the detection of novel species in the future.IMPORTANCEIn recent years, Campylobacter has gained recognition as the leading cause of bacterial gastroenteritis worldwide, leading to a substantial rise in the collection of genomic data of the Campylobacter genus in public databases. Currently, a standardized Campylobacter species boundary at the genomic level is absent, leading to challenges in detecting emerging pathogens and defining putative novel species within this genus. We used a comprehensive representation of genomes of the Campylobacter genus to construct a core genome phylogenetic tree. Furthermore, we found an average nucleotide identity (ANI) of 94.2% as the optimal cutoff to define the Campylobacter species. Using this cutoff, we identified 60 ANI genomic species which provided a standardized species definition and nomenclature. Importantly, we have developed Campylobacter Genomic Species typer (CampyGStyper), which can robustly and accurately assign these ANI genomic species to Campylobacter genomes, thereby aiding pathogen surveillance and facilitating evolutionary and epidemiological studies of existing and emerging pathogens in the genus Campylobacter.

Emergence of Poultry-Associated Human Salmonella enterica Serovar Abortusovis Infections, New South Wales, Australia.

Salmonella enterica serovar Abortusovis is a ovine-adapted pathogen that causes spontaneous abortion. Salmonella Abortusovis was reported in poultry in 2009 and has since been reported in human infections in New South Wales, Australia. Phylogenomic analysis revealed a clade of 51 closely related isolates from Australia originating in 2004. That clade was genetically distinct from ovine-associated isolates. The clade was widespread in New South Wales poultry production facilities but was only responsible for sporadic human infections. Some known virulence factors associated with human infections were only found in the poultry-associated clade, some of which were acquired through prophages and plasmids. Furthermore, the ovine-associated clade showed signs of genome decay, but the poultry-associated clade did not. Those genomic changes most likely led to differences in host range and disease type. Surveillance using the newly identified genetic markers will be vital for tracking Salmonella Abortusovis transmission in animals and to humans and preventing future outbreaks.

Genomic Characteristion of Opportunistic Pathogen Kluyvera Reveals a Novel CTX-M Subgroup.

A rising incidence of clinical infections has been caused by Kluyvera, a significant opportunistic pathogen. Meanwhile, Kluyvera acts as an important reservoir of blaCTX-Ms, which are the dominant genes of class A extended-spectrum ß-lactamases (ESBLs). In this work, 60 strains of Kluyvera were subjected to phylogenetic relationship reconstruction, antimicrobial susceptibility testing, and antibiotic resistance genes prediction. All mature blaCTX-Ms were gathered to perform subgroup reclassification. The findings demonstrate that Kluyvera has a large gene pool with significant genetic flexibility. Notably, 25% of strains showed simultaneous detection of ESBLs and carbapenem resistance genes. The genotypes of fourteen novel blaCTX-Ms were identified. A new subgroup classification approach for blaCTX-Ms was defined by using 20 amino acid site variants, which could split blaCTX-Ms into 10 subgroups. The results of the subgroup division were consistent with the phylogenetic clustering. More significantly, we proposed a novel blaCTX-M subgroup, KLUS, that is chromosomally encoded in K. sichuanensis and the new species put forward in this study, showing amino acid differences from the currently known sequences. Cloning and transformation tests demonstrated that the recipient bacteria had a robust phenotype of cefotaxime resistance. Closely related Kluyvera species had blaCTX-Ms in the same subgroup. Our research lays the groundwork for a deeper comprehension of Kluyvera and emphasizes how important a blaCTX-M reservoir it is. We provide an update on blaCTX-M subgroups reclassification from the aspects of phylogenetic relationship, amino acid differences, and the new subgroup KLUS, which needs to be strengthen monitored due to its strong resistance phenotype to cefotaxime.

Emergence and genomic insights of non-pandemic O1 Vibrio cholerae in Zhejiang, China.

IMPORTANCE: It is well recognized that only Vibrio cholerae O1 causes cholera pandemics. However, not all O1 strains cause pandemic-level disease. In this study, we analyzed non-pandemic O1 V. cholerae isolates from the 1960s to the 1990s from China and found that they fell into three lineages, one of which shared the most recent common ancestor with pandemic O1 strains. Each of these non-pandemic O1 lineages has unique properties that contribute to their capacity to cause cholera. The findings of this study enhanced our understanding of the emergence and evolution of both pandemic and non-pandemic O1 V. cholerae.

Emerging Aeromonas enteric infections: their association with inflammatory bowel disease and novel pathogenic mechanisms.

Aeromonas species are emerging human enteric pathogens. This study examines the isolation of Aeromonas and other enteric bacterial pathogens from patients with and without inflammatory bowel disease (IBD). This study also investigates the intestinal epithelial pathogenic mechanisms of Aeromonas veronii. The isolation rates of seven enteric bacterial pathogens from 2,279 patients with IBD and 373,276 non-IBD patients were compared. An A. veronii strain (AS1) isolated from intestinal biopsies of a patient with IBD was used for pathogenic mechanism investigation, and Escherichia coli K12 was used as a bacterial control. HT-29 cells were used as a model of human intestinal epithelium. A significantly higher isolation of Aeromonas species was found in patients with IBD as compared to non-IBD patients (P = 0.0001, odds ratio = 2.11). A. veronii upregulated 177 inflammatory genes and downregulated 52 protein-coding genes affecting chromatin assembly, multiple small nuclear RNAs, multiple nucleolar RNAs, and 55 cytoplasmic tRNAs in HT-29 cells. These downregulation effects were unique to A. veronii and not observed in HT-29 cells infected with E. coli K12. A. veronii induced intestinal epithelial apoptosis involving the intrinsic pathway. A. veronii caused epithelial microvilli shortening and damage and epithelial production of IL-8. In conclusion, this study for the first time reports the association between IBD and Aeromonas enteric infection detected by bacterial cultivation. This study also reports that A. veronii damages intestinal epithelial cells via multiple mechanisms, of which the downregulating cytoplasmic tRNA, small nuclear RNA, and small nucleolar RNA are novel bacterial pathogenic mechanisms. IMPORTANCE This study for the first time reports the association between inflammatory bowel disease (IBD) and Aeromonas enteric infection detected by bacterial pathogen cultivation, highlighting the need of clinical and public health attention. The finding that patients with IBD are more susceptible to Aeromonas enteric infection suggests that detection of Aeromonas enteric infection should be routinely performed for the diagnosis and treatment of IBD. This study also reports novel bacterial pathogenic mechanisms employed by Aeromonas veronii. Through comparative transcriptomic analysis and other techniques, this study revealed the pathogenic mechanisms by which A. veronii causes damage to intestinal epithelial cells. Among the various pathogenic mechanisms identified, the downregulating tRNA, small nuclear and nucleolar RNAs in human intestinal epithelial cells are novel bacterial pathogenic mechanisms.

Integrating proteomic data with metabolic modeling provides insight into key pathways of Bordetella pertussis biofilms.

Pertussis, commonly known as whooping cough is a severe respiratory disease caused by the bacterium, Bordetella pertussis. Despite widespread vaccination, pertussis resurgence has been observed globally. The development of the current acellular vaccine (ACV) has been based on planktonic studies. However, recent studies have shown that B. pertussis readily forms biofilms. A better understanding of B. pertussis biofilms is important for developing novel vaccines that can target all aspects of B. pertussis infection. This study compared the proteomic expression of biofilm and planktonic B. pertussis cells to identify key changes between the conditions. Major differences were identified in virulence factors including an upregulation of toxins (adenylate cyclase toxin and dermonecrotic toxin) and downregulation of pertactin and type III secretion system proteins in biofilm cells. To further dissect metabolic pathways that are altered during the biofilm lifestyle, the proteomic data was then incorporated into a genome scale metabolic model using the Integrative Metabolic Analysis Tool (iMAT). The generated models predicted that planktonic cells utilised the glyoxylate shunt while biofilm cells completed the full tricarboxylic acid cycle. Differences in processing aspartate, arginine and alanine were identified as well as unique export of valine out of biofilm cells which may have a role in inter-bacterial communication and regulation. Finally, increased polyhydroxybutyrate accumulation and superoxide dismutase activity in biofilm cells may contribute to increased persistence during infection. Taken together, this study modeled major proteomic and metabolic changes that occur in biofilm cells which helps lay the groundwork for further understanding B. pertussis pathogenesis.

Vibrio metschnikovii as an emergent pathogen: analyses of phylogeny and O-antigen and identification of possible virulence characteristics.

Vibrio metschnikovii is an emergent pathogen that causes human infections which may be fatal. However, the phylogenetic characteristics and pathogenicity determinants of V. metschnikovii are poorly understood. Here, the whole-genome features of 103 V. metschnikovii strains isolated from different sources are described. On phylogenetic analysis V. metschnikovii populations could be divided into two major lineages, defined as lineage 1 (L1) and 2 (L2), of which L1 was more likely to be associated with human activity. Meanwhile, we defined 29 V. metschnikovii O-genotypes (VMOg, named VMOg1-VMOg29) by analysis of the O-antigen biosynthesis gene clusters (O-AGCs). Most VMOgs (VMOg1 to VMOg28) were assembled by the Wzx/Wzy pathway, while only VMOg29 used the ABC transporter pathway. Based on the sequence variation of the wzx and wzt genes, an in silico O-genotyping system for V. metschnikovii was developed. Furthermore, nineteen virulence-associated factors involving 161 genes were identified within the V. metschnikovii genomes, including genes encoding motility, adherence, toxins, and secretion systems. In particular, V. metschnikovii was found to promote a high level of cytotoxicity through the synergistic action of the lateral flagella and T6SS. The lateral flagellar-associated flhA gene played an important role in the adhesion and colonization of V. metschnikovii during the early stages of infection. Overall, this study provides an enhanced understanding of the genomic evolution, O-AGCs diversity, and potential pathogenic features of V. metschnikovii.

Genomic epidemiology and multilevel genome typing of Bordetella pertussis.

Bordetella pertussis causes pertussis (or whooping cough), a severe respiratory infectious disease in infants, although it can be prevented by whole cell and acellular vaccines. The recent pertussis resurgence in industrialised countries is partly attributed to pathogen adaptation to vaccines, while emergence of antimicrobial resistance, specifically to macrolides in China, has become a concern. Surveillance of current circulating and emerging strains is therefore vital to understand the risks they pose to public health. Although the use of genomics-based typing is increasing a genomic nomenclature for this pathogen has not been well established. Here, we implemented the multilevel genome typing (MGT) system for B. pertussis with five levels of resolution, which provide targeted typing of relevant lineages and discrimination of closely related strains at the finest scale. The lower resolution levels (MGT2 and MGT3) describe the distribution of major vaccine antigen alleles including ptxP, fim3, fhaB and prn, as well as temporal and spatial trends within the B. pertussis global population. Mid-resolution levels (MGT3 and MGT4) enable typing of antibiotic-resistant lineages and Prn deficient lineages within the ptxP3 clade. The high-resolution level (MGT5) can capture finer-scale epidemiology such as outbreaks and local transmission events, with comparable resolution to existing genomic methods of strain-relatedness assessment. The scheme offers stable MGT-type assignments aiding harmonisation of typing and communication between laboratories. The scheme is available at https://mgtdb.unsw.edu.au/pertussis, is regularly updated from global data repositories and accepts public submissions. The MGT scheme provides a comprehensive, robust, and scalable system for global surveillance of B. pertussis.

Bacteroides vulgatus Ameliorates Lipid Metabolic Disorders and Modulates Gut Microbial Composition in Hyperlipidemic Rats.

Hyperlipidemia is a risk factor and key indicator for cardiovascular diseases, and the gut microbiota is highly associated with hyperlipidemia. Bacteroides vulgatus is a prevalent mutualist across human populations and confers multiple health benefits such as immunoregulation, antiobesity, and coronary artery disease intervention. However, its role in antihyperlipidemia has not been systematically characterized. This study sought to identify the effect of B. vulgatus Bv46 on hyperlipidemia. Hyperlipidemic rats were modeled by feeding them a high-fat diet for 6 weeks. The effect of B. vulgatus Bv46 supplementation was evaluated by measuring anthropometric parameters, lipid and inflammation markers, and the liver pathology. Multi-omics was used to explore the underlying mechanisms. The ability of B. vulgatus Bv46 to produce bile salt hydrolase was confirmed by gene annotation and in vitro experiments. Oral administration of B. vulgatus Bv46 in hyperlipidemic rats significantly reduced the body weight gain, food efficiency, and liver index, improved the serum lipid profile, lowered the levels of serum inflammatory cytokines, promoted the loss of fecal bile acids (BAs), and extended the fecal pool of short-chain fatty acids (SCFAs), especially propionate and butyrate. B. vulgatus Bv46 induced compositional shifts of the gut microbial community of hyperlipidemic rats, characterized by a lower ratio of Firmicutes to Bacteroidetes with an increase of genera Bacteroides and Parabacteroides. After intervention, serum metabolite profiling exhibited an adaptation in amino acids and glycerophospholipid metabolism. Transcriptomics further detected altered biological processes, including primary bile acid biosynthesis and fatty acid metabolic process. Taken together, the findings suggest that B. vulgatus Bv46 could be a promising candidate for interventions against hyperlipidemia. IMPORTANCE As a core microbe of the human gut ecosystem, Bacteroides vulgatus has been linked to multiple aspects of metabolic disorders in a collection of associative studies, which, while indicative, warrants more direct experimental evidence to verify. In this study, we experimentally demonstrated that oral administration of B. vulgatus Bv46 ameliorated the serum lipid profile and systemic inflammation of high-fat diet-induced hyperlipidemic rats in a microbiome-regulated manner, which appears to be associated with changes of bile acid metabolism, short-chain fatty acid biosynthesis, and serum metabolomic profile. This finding supports the causal contribution of B. vulgatus in host metabolism and helps to form the basis of novel therapies for the treatment of hyperlipidemia.

Genomic Epidemiology and Multilevel Genome Typing of Australian Salmonella enterica Serovar Enteritidis.

Salmonella enterica serovar Enteritidis is one of the leading causes of salmonellosis in Australia. In this study, a total of 568 S. Enteritidis isolates from two Australian states across two consecutive years were analyzed and compared to international strains, using the S. Enteritidis multilevel genome typing (MGT) database, which contained 40,390 publicly available genomes from 99 countries. The Australian S. Enteritidis isolates were divided into three phylogenetic clades (A, B, and C). Clades A and C represented 16.4% and 3.5% of the total isolates, respectively, and were of local origin. Clade B accounted for 80.1% of the isolates which belonged to seven previously defined lineages but was dominated by the global epidemic lineage. At the MGT5 level, three out of five top sequence types (STs) in Australia were also top STs in Asia, suggesting that a fair proportion of Australian S. Enteritidis cases may be epidemiologically linked with Asian strains. In 2018, a large egg-associated local outbreak was caused by a recently defined clade B lineage prevalent in Europe and was closely related, but not directly linked, to three European isolates. Additionally, over half (54.8%) of predicted multidrug resistance (MDR) isolates belonged to 10 MDR-associated MGT-STs, which were also frequent in Asian S. Enteritidis . Overall, this study investigated the genomic epidemiology of S. Enteritidis in Australia, including the first large local outbreak, using MGT. The open MGT platform enables a standardized and sharable nomenclature that can be effectively applied to public health for unified surveillance of S. Enteritidis nationally and globally. IMPORTANCE Salmonella enterica serovar Enteritidis is a leading cause of foodborne infections. We previously developed a genomic typing database (MGTdb) for S. Enteritidis to facilitate global surveillance of this pathogen. In this study, we examined the genomic features of Australian S. Enteritidis using the MGTdb and found that Australian S. Enteritidis is mainly epidemiologically linked with Asian strains (especially strains carrying antimicrobial resistance genes), followed by European strains. The first large-scale egg-associated local outbreak in Australia was caused by a recently defined lineage prevalent in Europe, and three European isolates in the MGTdb were closely related but not directly linked to this outbreak. In summary, the S. Enteritidis MGTdb open platform is shown to be a potentially powerful tool for national and global public health surveillance of this pathogen.

Bacteroides vulgatus attenuates experimental mice colitis through modulating gut microbiota and immune responses.

Introduction: Bacteroides vulgatus is one of the predominant Bacteroides species in the human gut and exerts a series of beneficial effects. The aim of this study was to investigate the protective role of B. vulgatus Bv46 in a dextran sodium sulfate (DSS) induced colitis mouse model. Methods: Female C57BL/6J mice were given 3% DSS in drinking water to induce colitis and simultaneously treated with B. vulgatus Bv46 by gavage for 7 days. Daily weight and disease activity index (DAI) of mice were recorded, and the colon length and histological changes were evaluated. The effects of B. vulgatus Bv46 on gut microbiota composition, fecal short chain fatty acids (SCFAs) concentration, transcriptome of colon, colonic cytokine level and cytokine secretion of RAW 264·7 macrophage cell line activated by the lipopolysaccharide (LPS) were assessed. Results and Discussion: B. vulgatus Bv46 significantly attenuated symptoms of DSS-induced colitis in mice, including reduced DAI, prevented colon shortening, and alleviated colon histopathological damage. B. vulgatus Bv46 modified the gut microbiota community of colitis mice and observably increased the abundance of Parabacteroides, Bacteroides, Anaerotignum and Alistipes at the genus level. In addition, B. vulgatus Bv46 treatment decreased the expression of colonic TNF-α, IL-1ß and IL-6 in DSS-induced mouse colitis in vivo, reduced the secretion of TNF-α, IL-1ß and IL-6 in macrophages stimulated by LPS in vitro, and downregulated the expression of Ccl19, Cd19, Cd22, Cd40 and Cxcr5 genes in mice colon, which mainly participate in the regulation of B cell responses. Furthermore, oral administration of B. vulgatus Bv46 notably increased the contents of fecal SCFAs, especially butyric acid and propionic acid, which may contribute to the anti-inflammatory effect of B. vulgatus Bv46. Supplementation with B. vulgatus Bv46 serves as a promising strategy for the prevention of colitis.

MGTdb: a web service and database for studying the global and local genomic epidemiology of bacterial pathogens.

Multilevel genome typing (MGT) enables the genomic characterization of bacterial isolates and the relationships among them. The MGT system describes an isolate using multiple multilocus sequence typing (MLST) schemes, referred to as levels. Thus, for a new isolate, sequence types (STs) assigned at multiple precisely defined levels can be used to type isolates at multiple resolutions. The MGT designation for isolates is stable, and the assignment is faster than the existing approaches. MGT's utility has been demonstrated in multiple species. This paper presents a publicly accessible web service called MGTdb, which enables the assignment of MGT STs to isolates, along with their storage, retrieval and analysis. The MGTdb web service enables upload of genome data as sequence reads or alleles, which are processed and assigned MGT identifiers. Additionally, any newly sequenced isolates deposited in the National Center for Biotechnology Information's Sequence Read Archive are also regularly retrieved (currently daily), processed, assigned MGT identifiers and made publicly available in MGTdb. Interactive visualization tools are presented to assist analysis, along with capabilities to download publicly available isolates and assignments for use with external software. MGTdb is currently available for Salmonella enterica serovars Typhimurium and Enteritidis and Vibrio cholerae. We demonstrate the usability of MGTdb through three case studies - to study the long-term national surveillance of S. Typhimurium, the local epidemiology and outbreaks of S. Typhimurium, and the global epidemiology of V. cholerae. Thus, MGTdb enables epidemiological and microbiological investigations at multiple levels of resolution for all publicly available isolates of these pathogens. Database URL: https://mgtdb.unsw.edu.au.

Genomic Epidemiology of Vibrio cholerae O139, Zhejiang Province, China, 1994-2018.

Cholera caused by Vibrio cholerae O139 was first reported in Bangladesh and India in 1992. To determine the genomic epidemiology and origins of O139 in China, we sequenced 104 O139 isolates collected from Zhejiang Province, China, during 1994-2018 and compared them with 57 O139 genomes from other countries in Asia. Most Zhejiang isolates fell into 3 clusters (C1-C3), which probably originated in India (C1) and Thailand (C2 and C3) during the early 1990s. Different clusters harbored different antimicrobial resistance genes and IncA/C plasmids. The integrative and conjugative elements carried by Zhejiang isolates were of a new type, differing from ICEVchInd4 and SXTMO10 by single-nucleotide polymorphisms and presence of genes. Quinolone resistance-conferring mutations S85L in parC and S83I in gyrA occurred in 71.2% of the Zhejiang isolates. The ctxB copy number differed among the 3 clusters. Our findings provided new insights for prevention and control of O139 cholera .

Evolutionary and genomic insights into the long-term colonization of Shigella flexneri in animals.

The enteroinvasive bacterium Shigella flexneri is known as a highly host-adapted human pathogen. There had been no known other reservoirs reported until recently. Here 34 isolates obtained from animals (yaks, dairy cows and beef cattle) from 2016 to 2017 and 268 human S. flexneri isolates from China were sequenced to determine the relationships between animal and human isolates and infer the evolutionary history of animal-associated S. flexneri. The 18 animal isolates (15 yak and 3 beef cattle isolates) in PG1 were separated into 4 lineages, and the 16 animal isolates (1 yak, 5 beef cattle and 10 dairy cow isolates) in PG3 were clustered in 8 lineages. The most recent human isolates from China belonged to PG3 whereas Chinese isolates from the 1950s-1960s belonged to PG1. PG1 S. flexneri may has been transmitted to the yaks during PG1 circulation in the human population in China and has remained in the yak population since, while PG3 S. flexneri in animals were likely recent transmissions from the human population. Increased stability of the large virulence plasmid and acquisition of abundant antimicrobial resistance determinants may have enabled PG3 to expand globally and replaced PG1 in China. Our study confirms that animals may act as a reservoir for S. flexneri. Genomic analysis revealed the evolutionary history of multiple S. flexneri lineages in animals and humans in China. However, further studies are required to determine the public health threat of S. flexneri from animals.

Genomic dissection of the microevolution of Australian epidemic Bordetella pertussis.

ABSTRACTWhooping cough (pertussis) is a highly contagious respiratory disease caused by the bacterium Bordetella pertussis. Despite high vaccine coverage, pertussis has re-emerged in many countries including Australia and caused two large epidemics in Australia since 2007. Here, we undertook a genomic and phylogeographic study of 385 Australian B. pertussis isolates collected from 2008 to 2017. The Australian B. pertussis population was found to be composed of mostly ptxP3 strains carrying different fim3 alleles, with ptxP3-fim3A genotype expanding far more than ptxP3-fim3B. Within the former, there were six co-circulating epidemic lineages (EL1 to EL6). The multiple ELs emerged, expanded, and then declined at different time points over the two epidemics. In population genetics terms, both hard and soft selective sweeps through vaccine selection pressures have determined the population dynamics of Australian B. pertussis. Relative risk estimation suggests that once a new B. pertussis lineage emerged, it was more likely to spread locally within the first 1.5 years. However, after 1.5 years, any new lineage was likely to expand to a wider region. Phylogenetic analysis revealed the expansion of ptxP3 strains was also associated with replacement of the type III secretion system allele bscI1 with bscI3. bscI3 is associated with decreased T3SS secretion and may allow B. pertussis to reduce immune recognition. This study advanced our understanding of the epidemic population structure and spatial and temporal dynamics of B. pertussis in a highly immunized population.

Dissecting Listeria monocytogenes Persistent Contamination in a Retail Market Using Whole-Genome Sequencing.

Listeria monocytogenes is a foodborne pathogen that can cause invasive disease with high mortality in immunocompromised individuals and can survive in a variety of food-associated environments for a long time. L. monocytogenes clonal complex (CC) 87 is composed of ST87 and three other STs and has been identified as the most common subgroup associated with both foods and human clinical infections in China. Therefore, the persistence of CC87 L. monocytogenes in food-associated environments poses a significant concern for food safety. In this study, 83 draft genomes of CC87 L. monocytogenes, including 60 newly sequenced genomes, were analyzed with all isolates from our previous surveillance in Zigong, Sichuang, China. Sixty-eight of the studied isolates were isolated from one retail market (M1 market), while the others were from seven other markets (M2-M8 markets) in the same city. Whole-genome multilocus sequence typing (wg-MLST) and the whole-genome single nucleotide polymorphism (wg-SNP) analysis were performed. Three persistent contamination routes were identified in the M1 market, caused by 2 clusters (A and B) and a wgST31 type. Cluster A isolates were associated with the persistent contamination in a raw meat stall (M1-S77), while Cluster B isolates caused a persistent contamination in aquatic foods stalls. Five wgST31 isolates caused persistent contamination in a single aquatic stall (M1-S65). A pLM1686-like plasmid was found in all Cluster A isolates. A novel plasmid, pLM1692, a truncated pLM1686 plasmid without the cadmium, and other heavy metal resistance genes were conserved in all wgST31 isolates. By comparing persistent and putative non-persistent isolates, four genes that were all located in the prophage comK might be associated with persistence. These findings enhanced our understanding of the underlying mechanisms of contamination and assist in formulating targeted strategies for the prevention and control of L. monocytogenes transmission from the food processing chain to humans. IMPORTANCE Contamination of food by Listeria monocytogenes at retail level leads to potential consumption of contaminated food with high risk of human infection. Our previous study found persistent contamination of CC87 L. monocytogenes from a retail market in China through pulsed-field gel electrophoresis and multilocus sequence typing. In this study, whole-genome sequencing was used to obtain the highest resolution inference of the source and reasons for persistent contamination; meat grinders and minced meat were the major reservoir of persistent contamination in meat stalls, whereas fishponds were the major reservoir in seafood stalls, with different L. monocytogenes isolates involved. These isolates carried different properties such as plasmids and prophages, which may have contributed to their ability to survive or adapt to the different environments. Our findings suggest that whole-genome sequencing will be an effective surveillance tool to detect persistent L. monocytogenes contamination in retail food markets and to design new control strategies to improve food safety.

Horizontal transfer of antibiotic resistance genes within the bacterial communities in aquacultural environment.

Very little is known about how microbiome interactions shape the horizontal transfer of antibiotic resistance genes in aquacultural environment. To this end, we first conducted 16S rRNA gene amplicon sequencing to monitor the dynamics of bacterial community compositions in one shrimp farm from 2019 to 2020. Next, co-occurrence analysis was then conducted to reveal the interactions network between Vibrio spp. and other species. Subsequently, 21 V. parahaemolyticus isolates and 15 related bacterial species were selected for whole-genome sequencing (WGS). The 16S rDNA amplicon sequencing results identified a remarkable increase of Vibrio and Providencia in September-2019 and a significant rise of Enterobacter and Shewanella in Septtember-2020. Co-occurrence analysis revealed that Vibrio spp. positively interacted with the above species, leading to the sequencing of their isolates to further understand the sharing of the resistant genomic islands (GIs). Subsequent pan-genomic analysis of V. parahaemolyticus genomes identified 278 horizontally transferred genes in 10 GIs, most of which were associated with antibiotic resistance, virulence, and fitness of metabolism. Most of the GIs have also been identified in Providencia, and Enterobacter, suggesting that exchange of genetic traits might occur in V. parahaemolyticus and other cooperative species in a specific niche. No genetic exchange was found between the species with negative relationships. The knowledge generated from this study would greatly improve our capacity to predict and mitigate the emergence of new resistant population and provide practical guidance on the microbial management during the aquacultural activities.